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Abstract 

Here we comment on the thermodynamic inconsistency problem and the reformulation 
of statistical mechanics of widely studied quasiparticle models of quark-gluon plasma. Their 
starting relation, the expression for pressure itself is a wrong choice and lead to thermo- 
dynamic inconsistency and the requirements of the reformulation of statistical mechanics. 
We propose a new approach to the problem using the standard statistical mechanics and is 
thermodynamically consistent. We also show that the other quasiparticle models may be 
obtained from our general formalism as a special case under certain restrictive condition. 
Further, as an example, we have applied our model to explain the nonideal behaviour of 
gluon plasma and obtained a remarkable good fit to the lattice results by adjusting just a 
single parameter. 
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1 Introduction 



The quasiparticle quark-gluon plasma (qQGP) is a phenomenological model, with few fitting pa- 
rameters, widely used to describe the nonideal behaviour of quark-gluon plasma (QGP). It was 
first proposed by Goloviznin and Satz [1] and then by Peshier et. al. [2] to explain the equation 
of state (EoS) of QGP from lattice gauge theory (LGT) simulation of quantum chromodynamics 
(QCD) at finite temperature [3]. The model, however, failed [5] to explain the more accurate, 
recent lattice data [4]. Further, Gorenstein and Yang [6] pointed out that the model is thermo- 
dynamically inconsistent. This thermodynamically inconsistency problem was remedied by them 
by introducing a temperature dependent vacuum energy and forced it to cancel the thermody- 
namically inconsistent term, which was named as the reformulated statistical mechanics. It is 
still not clear what is the physics or origin of this constraint which was called as thermodynamic 
consistency check in Ref. [6, 7, 8, 9] Here we show that the whole exercise is unnecessary and 
following the standard statistical mechanics (SM), we propose a new qQGP model which contains 
a single phenomenological parameter. Our model is thermodynamically consistent and explains 
lattice data very well. 

2 Our model of qQGP: 

Let us start with the work of Peshier et. al. [2] on gluon plasma. All thermodynamic quantities 
were derived from the pressure, P, which was assumed as 
PV 00 

_ = -5>(i- e -* fc)> (i) 

where the right hand side is the logarithm of the grand partition function, Qq(T), and e& is the 
single particle energy of quasi-gluon, i.e, gluon with temperature dependent mass, given by, 

e k = y/k 2 + m?(T) , 

where k is momentum and m is mass. (3 is defined as 1/T. The expression for pressure is similar 
to that of ideal gas with temperature dependent mass and hence let us denote it as P^. We 
want to stress that this assumption itself is the root cause of thermodynamic inconsistency and 
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hence the reformulation of SM by Gorenstein and Yang [6] . Generally, in grand canonical ensemble 
(GCE), energy (E r ) and number of particles (N s ) fluctuate, but temperature (T) and the chemical 
potential (//) are fixed. Hence, the average energy (U) and average number of particles (N) are 
defined and may be related to the grand partition function or q-potential, 

oo 

q = ln^G = HT, e ~ PEr ~ aNs ) =TE Tze~^) , (2) 

s,r k=0 

where =F for bosons and fermions and z = e^ T = e~ a is called fugacity. The average energy U is 
defined as, 

Qg df3 Y 1 =F z e 



Us<Er>=i ^ = ( 3) 



Note that the partial differentiation with repect to f3 above is just a mathematical method to 
express U in terms of sum over single particle energy levels, e^, making use of Eq. (2). While 
differentiating, indirect dependence of j3 = 1/T in the fugacity, z, and mass, m(T), must be 
ignored by definition. Otherwise, we will not get back < E r >. Similarly, the average density N 
is defined as, 

These (Eqs. (3), (4)) are the standard relations [10] of U and N to the partition function, which 
is valid even for quasiparticle with (T,/x) dependent mass by the definition of averages. Here, for 
gluon plasma, we have // = or z — 1. Next, pressure may be obtained by two methods. In 
method-I, one starts from U and using thermodynamic relation, 

^l = T d ^-P, (5) 
V dT K J 

and on integration, one gets pressure which is the procedure that we follow here. In method II, 
again following the standard text books on SM [10], we can relate P to q-potential as follows. 
The variation in q-potential due to variations in it's dependence, namely T, ji and volume V, 
specifying the macro-state of GCE system, is, 



(6) 
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Now, when compared with the text books results, we have an extra term coming from 5E r due to 
temperature dependent mass and then using the definition of averages, we get, 



Therefore we see that P is not just equal to P id , but there is an extra term. This extra term 
ensure thermodynamic consistency of the relation as follows. From above P, on differentiating 
with respect to T for a system with /i — or z — 1, 

&T~T + f~V < df >+ V < ~dT > ^ 
where the last two terms exactly cancels (following the procedure used in deriving Eq, (3)) and 
hence the thermodynamic relation, Eq. (5), is obeyed as expected. 

Further more, this P is also consistent with the P obtained from U through thermodynamic 
relations which may be shown as follows. The Eq. (7) may be simplified by evaluating < ^ > 
and taking continuum limit (V — > oo), for a system with /i — 0, as, 

which on simplification, reduces to, 
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(10) 



where gf is the internal degrees of freedom and K ± , K 2 are modified Bessel functions. Using the 
recurrence relations of Bessel functions and on integration by parts, above equation may be further 
simplified to get, 

P Pn fl 3 

r=^~ dps, (11) 

1 i J/3 

where e is the energy density and Pq is the pressure at some temperature To or Pq. This equation 
is nothing but the thermodynamic relation, Eq. (5). Therefore, Gorenstein and Yang's starting 
argument that above two methods give different m(T) does not exist now by using our derived 
expression for P, instead of the assumption [2, 6]. 



3 Question of vacuum energy B(T) : 

After the reformulation of SM by Gorenstein and Yang, almost all study in qQGP is based on 
the thermodynamic consistency relation, related to vacuum energy B(T). Different authors call 
and interpret B(T) in a different way, like vacuum energy, background field or bag pressure. But, 
by definition of quasiparticle, whole thermal energy is used to excite these quasiparticles. So 
quasiparticles are excitations above the ground state or vacuum state which may not depend on 
temperature or chemical potential. This is our assumption. As noted earlier, we also don't have 
any thermodynamic inconsistency in our model. 

In fact, when we redo our derivation of pressure with vacuum or zero point energy in single 
particle energy, like in Ref. [6], Eq. (9) is modified as, 

/ r T dr 

P = P ld -B{T)+T [ - 

\Jt t 

and the energy density, 

e = e td + B(T). (13) 

where is the expression for e similar to ideal gas. Again it is easy to show that above P 
and e obey thermodynamic relation Eq. (5). The thermodynamic consistency relation [6], used 
in other qQGP models, is nothing but a restrictive condition that the terms inside the square 
bracket in Eq. (12) is zero. At present it is not clear what is the physical origin of this constraint. 
Note that without this constraint, so called thermodynamic consistency relation, our system is 
thermodynamically consistent even with the zero-point energy contribution, B(T). One may model 
B(T) based on other effects of strongly interacting QCD system, like hadronic states, resonances 
and may be relevent at the transition point. In our study of gluon plasma here, we neglect all 
these effects and consider only the thermal properties of gluons. Hence we take B(T) = and 
we get a very good fit to lattice results except at very close to the transition temperature, i.e, for 
T < 1.2T C . 



gf dm 



— / dk 
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dB 



e k (e^ T 1) dr 



(12) 
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4 EoS of gluon plasma: 



As an example, let us apply our model to gluon plasma which is a QCD plasma without quarks. 
We first calculate the energy density, expressed in terms of e(T) = e/e s , and then obtain P from 
thermodynamic relation, Eq. (5). So we have, from Eq. (3) after some algebra, 



1 r oo -, 



(14) 



where e s is the Stefan-Boltzman gas limit of QGP, m g is the temperature dependent mass and K ± , 



K 2 are modified Bessel functions. The results are plotted in Fig. 1, for two different mass terms, 
m? g (T) = u) 2 p = g 2 (T)T 2 /3 (our model) and m 2 g (T) = g 2 (T)T 2 /2 (other qQGP models). g 2 (T) is 
related to the two- loop order running coupling constant, given by, 

„ m _ ^ L 3(153 - lQn,) ln(2 ln(T/A r )) \ 

A) (33 - 2n f ) m(T/A T ) \ (33 - 2n 7 ) 2 ln(T/A T ) J' 1 ' 

where A^ is a parameter related to QCD scale parameter. This choice of a s (T) is an approximate 
expression of the running coupling constant used in lattice simulations [4] . Then the pressure is 
obtained from the thermodynamic relation Eq. (5) or Eq. (11). Since we have only one parameter 
to adjust, we don't get good fit for the generally used second choice of quasi-gluon mass. The 
best fitted parameter is A T /T C = 0.3. But a remarkably good fit may be obtained for our choice 
of gluon mass which is motivated from the fact that the quasi-gluons are the thermal excitations 
of plasma waves with mass equal to the plasma frequency [11]. The value of the fitted parameter 
is A T /T C = 0.65. 

Let us now compare our results with the results from other qQGP models, for example Ref. 
[7], where B(T) is not zero, but is determined by thermodynamic consistency relation. From the 
Fig. 2, we see that only at large T/T c both the results almost match, but differ near to T/T c = 1. 
We used the same ct s (T) with A T /T C = .65 for both the cases. Further, results from our model 
with B(T) = fits well the lattice data. A very good fit to lattice data was also obtained in 
Ref. [7], but with a different expression for a s (T), having two free parameters, and an additional 
parameter related to degrees of freedom. 
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5 Conclusions: 



Here we have pointed out the basic reason for the thermodynamic inconsistency of the extensively 
studied quasiparticle QGP models [2] and it's consequence of the reformulation of statistical me- 
chanics [6]. To revise it, we have proposed a new qQGP model which follows from the standard SM 
and has no thermodynamic inconsistency. When we extend our formalism to a system with tem- 
perature dependent vacuum energy, again, we get a thermodynamically consistent general model 
and we obtained other widely studied qQGP models as a special case of our model under certain 
restrictive condition which was called as thermodynamic consistency relation in Ref. [6, 7, 8, 9]. 
As an example, we studied the gluon plasma using our model. A remarkable good fit to the 
LGT data was obtained by adjusting just one parameter and without the temperature dependent 
vacuum energy B(T). Whereas we know that the other qQGP models has 3 or more parameters. 
Further extension of our model to flavored QGP without and with masses, and also without and 
with chemical potential, fit remarkably well the lattice results and were reported in [12, 13]. 
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Figure 1: Plots of P/T 4 (lower set of graphs) and e/T 4 (upper set of graphs) as a function of T/T c 
from our model and lattice results (symbols) [4] for gluon plasma with two different models for 
mass, m 2 (T) = g 2 T 2 /3 (dashed line) and m 2 {T) = g 2 T 2 /2 (dashed-dotted line). 

Figure 2: Plots of P/T 4 (lower set of graphs) and e/T 4 (upper set of graphs) as a function of T/T c 
from our model and lattice results (symbols) [4] for gluon plasma with B(T) = (dashed line) 
and using thermodynamic consistency relation (dashed-dotted line). 
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